Statistics of branched flow in a weak correlated random potential 
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Recent images of electron flow through a two-dimensional electron gas (2DEG) device show 
branching behavior that is reproduced in numerical simulations of motion in a correlated random 
potential [3]. We show how such branching naturally arises from caustics in the classical flow and 
find a simple scaling behavior of the branching under variation of the random potential strength. 
Analytic results describing the statistical properties of the branching are confirmed by classical and 
quantum numerical tests. 
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During the past decade, major advances in the fabri- 
cation and study of mesoscopic (or nanoscale) structures 
have led to the observation of many novel phenomena, 
and a fruitful interaction has resulted with simultane- 
ous theoretical progress in the areas of random matrix 
theory, quantum chaos, and disordered systems. Exam- 
ples of phenomena that have received much attention 
include universal conductance fluctuations and weak lo- 
calization in open quantum dots, and the statistics of 
conductance peak spacings and heights in the Coulomb 
blockade regime [|J. 

More recently, scanning probe microscope technology 
has allowed for the imaging of current flowing in a 2DEG 
device ||. After the electron current passes through a 
narrow quantum point contact (QPC), the flow shows 
a striking branch- like behavior [pi. The qualitative fea- 
tures of the observed branching pattern are well repro- 
duced in numerical simulations of quantum wave evolu- 
tion through a correlated random potential, when the 
rms potential height and the correlation length of the 
potential are fixed at their experimentally measured val- 
ues. Perhaps more surprisingly, the same branches are 
observed in the corresponding classical simulation, indi- 
cating that the experimental phenomenon should have 
a classical explanation. However, the simulations show 
that the branches do not correspond to valleys in the ran- 
dom potential. Instead, it was suggested that branches 
may arise from caustics in the classical flow. 

We emphasize that the problem of wave flow through 
a correlated random potential is of general interest and 
not restricted to the mesoscopic context. For example, 
similar effects of caustics have been discussed recently in 
the context of long-range sound propagation through the 
ocean Q. We also note that the influence of caustics 
on waves in random media has been studied as far back 
as 1977 by Berry who analyzed the moments of the 
intensity distribution. The related problem of the effect 
of orbit bifurcations on chaotic wave functions is also a 
research area of continuing interest || . Accordingly, our 
goal here is not only to make contact with one particu- 
lar experimental observation but also to make progress 
more generally in exploring this regime of classical and 
quantum chaotic behavior. 

In the present work, we extend the previous qualita- 
tive understanding of branching behavior, arriving at an- 
alytical results to be compared with numerics and ulti- 



mately with experiment. We focus on propagation in a 
two-dimensional potential, though the methods apply as 
well to the three-dimensional case (where branches are re- 
placed by caustic surfaces) and to related problems such 
as light propagation through a medium with varying re- 
fractive index. We concentrate throughout on generic 
patterns and not on rare events that would occur only in 
a few realizations of the random ensemble. 

Our starting point is a Gaussian-distributed random 
potential V(x,y) with a Gaussian correlation function, 



V{r) = 



V(f)V(r') = v%t 



(1) 



Without loss of generality, we will use units in which the 
correlation scale £ = 1 ; similarly the mass and the initial 
momentum, taken to be in the x— direction, are set to 
unity. To prevent backscattering and localization on a 
scale of order £, we must require the random potential 
to be weak relative to the initial kinetic energy, vo <C 
Ef = 1/2. In the experiment j|, vq/Ef w 0.08, and 
all our subsequent analytical expressions are computed 
to leading order in vq; however, the results turn out to 
be qualitatively and even in some respects quantitatively 
valid for vq/Ef = 0.3 and higher. In this approximation, 
then, the motion is unidirectional: p x (t) — p x (0) = 1, 
x(t) = t, and 

dy(t)/dt = p y (t) , d Py (t)/dt = -dV(t, y)/dy . (2) 

In other words, the two-dimensional dynamics reduces to 
one-dimensional evolution in the transverse dimension y, 
under the influence of an effectively time- varying random 
potential 0. For numerical purposes, the above classi- 
cal or quantum y-cvolution may be discretized on a scale 
At = Ax <C 1. The initial condition is zero transverse 
momentum; classically it is the manifold p y = in the 
y — p y phase space. 

Stretching and folding of this initial manifold under 
classical evolution in t (or x) produces simple caustics, 
or folds, where 5p y (t)/Sy(Q) — oo and the classical in- 
tensity in position space diverges |^| . By calculating how 
long it takes for a trajectory to travel a distance of one 
correlation length in the transverse y— direction, we may 
easily check that this stretching and folding occurs in the 
classical dynamics on the characteristic time scale 
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(3) 



i.e., when the potential is weak, a trajectory must pass 
over many correlation lengths of the potential in the lon- 
gitudinal (x) direction before caustics arise. In the fol- 
lowing, wc eliminate the explicit dependence of branching 
statistics on the random potential strength vq by express- 
ing all results in terms of the scaled time t/to. 

For t/to > 1) the transverse stretching factor defined 
as s = Sy(t) / Sy(0) develops a log-normal distribution, 



In P(s) 



(In s — at/tg) 2 



(4) 



where the exponent a is a dimensionless Lyapunov expo- 
nent, In s — at /to, while (3 characterizes the inhomogene- 
ity of the stretching, (51ns) 2 = (f3/2)t/t . The average 
amount of stretching in the y— direction of an infinitesi- 
mal piece of the initial manifold grows exponentially with 
time as 



6y(t)/5y(0)~e (a+l3/4)t/t ° 



(5) 



Since folds in the classical manifold typically develop 
when 8y ~ £ = 1, the number of caustics also grows 
exponentially with time as 



(6) 



where W is the length in the y— direction of the ini- 
tial manifold, i.e., the transverse width of the system in 
units of the correlation length. The scaling (^) has been 
checked numerically using several values of the random 
potential strength vq. 

Because the classical density of trajectories diverges 
at a caustic, we must perform smoothing of intensity on 
some scale b <C 1 in the transverse y— direction in or- 
der to obtain a well-defined branch height. In quantum 
mechanics, such smoothing is of course taken care of au- 
tomatically by the uncertainty principle, since caustics 
cannot be resolved when the phase space area enclosed 
by the fold is below U, and thus 6 c ff ~ fi 2 ^ 3 (see below). 
We also adopt the convention that the initial intensity 
(classical density of trajectories or quantum wave func- 
tion density) is normalized to unity. 

When a caustic forms at time t[ orm and transverse po- 
sition y, the classical 6— smoothed intensity i&(if rm, y) of 
the resulting branch scales as 



h 



-1/2 



s(0,tf OIm )]~ 



(J) 



where s(0, tf OTm ) is of course the stretching factor between 
time and time tf olm of the piece of the manifold in which 
the fold occurs. Away from caustics the intensity scales 
simply as ~ [s(0, i)] _1 . The extra factor of fo -1 / 2 in 
Eq. (Q) ensures that for small b (corresponding to the 
semiclassical or ray limit of the original quantum or wave 
problem) , caustics always dominate the tail of the inten- 
sity distribution, as we have confirmed numerically. 



Using the typical behavior of the stretching factor s 
(Eq. (||)), we now immediately see from Eq. (Q) that for 
or 1 < tfarmAo < |a _1 ln6 _1 , every caustic leads to a 
visible branch, with intensity 1. Once a given 

branch forms, its intensity immediately begins to de- 
cay (though at a slower rate) due to further stretching: 
h(t,y) ~ 6- 1/2 [s(0,t form )]- 1 [s(*form,t)]" 1/2 , eventually 
vanishing into the background on a logarithmic time scale 
t/to ~ a ln& ^> 1. The number of visible branches 
reaches its peak value of jV branch ~ Wb~ x l 2 at the time 
scale t/to ~ ^a -1 ln& _1 . 

At longer times, t/to > \ 0L ~ 1 m b^ 1 , it is no longer true 
that a typical new caustic results in a visible branch, since 
the resulting smoothed intensity is generally below the 
level (It, w 1) of the random background. Thus, at these 
longer times branches appear only at caustics that are 
formed from pieces of the manifold that have experienced 
anomalously little stretching compared with the average 
behavior. The number of such visible branches may be 
easily estimated by first calculating the number of caus- 
tics coming from those regions of the manifold that have 
stretched by no more than some factor so <C exp(at / to) ■ 
This number is given by the product of the manifold frac- 
tion that has stretched by less than sq and the number 
of caustics per unit length resulting from that stretching: 
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(8) 



Inserting the maximal stretching factor so = 
e {const /2)j ) -i/2 tnat w ju s till allow a newly-formed 

caustic to produce a visible branch, we obtain for 
t/to > \oT x Info" 1 
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(9) 



Thus, after a rapid initial climb, the number of branches 
falls off exponentially with scaled time t/to, at the rate 
a 2 1(3. The inherent ambiguity in defining how high a 
branch must be above the background to be visible can 
be absorbed into a multiplicative constant in JV branch and 
does not affect the falloff rate. 
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FIG. 1. Number of visible branches as a function of dimen- 
sionless scaled time t/to, where the characteristic time scale to 
depends on the strength wo of the random potential, Eq. (|^). 
A visible branch is defined as an intensity enhancement of at 
least a factor of 10 above the background. The two dashed 
curves correspond to different values of the smoothing length 
b. The solid line represents the predicted exponential falloff 
at long times (Eq. (|)), with slope a 2 / (5 = 0.30 (a = 0.65, 
P = 1.42). 

Numerical results are presented in Fig. |l|. In this fig- 
ure and throughout, a single realization of the random 
ensemble is used for each data curve; no ensemble aver- 
aging is necessary to obtain convincing statistics. 

In the experimental electron flow images ||, the few 
longest and most intense branches are the most visually 
striking. Let us therefore focus on the maximum branch 
intensity J™ ax as a function of scaled time t/t . This 
from Eq. (Q) is simply fe -1 / 2 divided by the minimum 
stretching factor, which is immediately found by setting 
JV^ = 1 in Eq. (||) and solving for lns . We obtain 



FIG. 2. Exponential falloff of maximum branch intensity 
7™ ax with scaled time t/to- The three dashed curves cor- 
respond to classical simulations for several values of vo and 
smoothing scale b, and the solid curve is a quantum calcu- 
lation. The power-law dependence of 7™ ax on b confirms 
that the peaks arise from simple caustics (folds). Finally, 
the straight line corresponds to the theoretical prediction 
of Eq. (llOj) , with analytically predicted exponent 7 = 0.17 
(a = 0.65, 13 = 1.42). 

Due to intensity fluctuations within a given branch, 
particularly in the quantum case, the fraction of space 
covered by branches, J branch 7 is in practice a more ro- 
bust measure than the branch number 7V branch (Eq. @). 
Defining J branch as the fraction of intensities that exceed 
some arbitrary cutoff J cutoff ; and combining the result of 
Eq. (0) with the fact that intensity falls off as the in- 
verse square root of the distance as we move away from 
a caustic, one straightforwardly obtains (for (3 > 2a) 

ln jbranch _ _ 2j ± _ 2 hi J Cutoff (11) 



ln/r^lnb- 1 -!^ 



This predicted exponential falloff of /™ ax with t/t is 
confirmed in Fig. g. We note the correct scaling behav- 
ior of the maximum intensity for different values of the 
potential strength vq and different smoothing lengths 6, 
as well as good agreement between the theoretically pre- 
dicted and numerically observed numerical exponent 7. 
The behavior 7™ ax ~ confirms that the maximum 

intensities indeed always come from caustics, which we 
have also confirmed by direct examination. 
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FIG. 3. Fraction of space J branch covered by branches taller 
than the arbitrary cutoff J cutoff = 10, as a function of scaled 
time t/to- Again, correct scaling is observed under changes 
of the potential strength vo . The slope of the solid line is the 
analytical prediction 27 = 0.34 of Eq. (p"l|). 

Eq. ([ll]) predicts that irrespectively of our definition of 
a branch (i.e., choice of J cutoff ) ; the spatial area covered 
by visible branches falls off exponentially with distance 
from the starting point, and the exponent is given by 27. 
This behavior is confirmed in Fig. |^ for several sets of 
parameters. We note in particular that to leading order 
the result does not depend on the smoothing scale b, and 
correspondingly on the ratio of wavelength to correlation 
length in the quantum or wave problem. Eq. ( |TT| ) may 
alternatively be interpreted as giving the cumulative dis- 
tribution of intensities at a fixed scaled time t/to: 

J™ dl' b P(Ii) ~ exp (-27^) h' 2 . (12) 
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i.e., P{h) ~ It. 3 - This power-law behavior remains 
valid until we reach the maximum intensity as given by 
Eq. @. The I~ 3 falloff means that tall branches do 
not affect the average intensity, but the mean squared 
intensity or inverse participation ratio (IPR) is logarith- 
mically divergent in the smoothing scale b. The predicted 
power-law intensity distribution at fixed time is tested 
numerically in Fig. ^. 
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FIG. 4. Cumulative probability distribution dl' b P(I' b ) 
of intensities I b at a fixed scaled time t/to = 2.2. The slope —2 
of the solid line is the theoretical prediction of Eq. (^|) . Note 
that the power-law behavior breaks down as we approach the 



maximum intensity from Eq. ([lo]) , 



e 4 5 for b = 0.0003 



and K 



for b = 0.00003. 



By explicitly tracking all caustic locations for the clas- 
sical dynamics in Fig. we have also checked that the 
fraction of large intensities lb that come from branches 
associated with caustics ranges from 98.7% for In/;, > 2 
to 99.8% for \nl b > 4, while only 3.8% of space overall is 
covered by such branches at time t/t^ — 2-2. This again 
provides confirmation that caustics dominate the tail of 
the intensity distribution. 

Our classical results for the branching statistics can be 
directly applied to the quantum context once we under- 
stand the smoothing scale b C R implied by the uncertainty 
principle in the quantum case. The parabolic structure 
of the classical manifold near a fold at [y* ,p y ) takes the 
form y - y* ~ [(p y - p* y )/p y yv ] 2 , where p y yp is the typi- 
cal value of the transverse momentum, p* yp ~ v^t 1 / 2 ~ 
(t/to) 1 / 2 . The area enclosed by the fold then scales as 

-y*fl 2 v 2 J 3 . Set- 



2/3 



(y-y*)(p y -p*y) ~ (y-y*) 3/2 p y yp ~ d/-yT"»o 

ting this area to h and dropping factors of order unity, 

2/3 4/9 

we finally obtain y — y* ~ h ' v , i.e., caustics in 
quantum mechanics do not get resolved below an effec- 
tive smoothing scale 



We may now simply transcribe our previously ob- 
tained classical results, substituting the effective quan- 
tum smoothing scale b c g for the explicit classical smooth- 
ing scale b. Data from typical quantum calculations ap- 
pear as solid curves in Figs. Q% 

In summary, we have made progress towards an an- 
alytic and quantitative understanding of branching be- 
havior for flow in a random correlated potential. The 
key conditions of applicability of these results are (1) 
a weak random potential vo/Ep <C 1, (2) a focus on 
distance scales of order the Lyapunov length, ^Lyap ~ 
1;(vq/Ef)~ 2 / 3 , where £ is the correlation scale of the po- 
tential, and (3) in the quantum case a wavelength no 
larger than £ to allow for a simple semiclassical treat- 
ment of the caustics. We mention in conclusion that 
the numerical values of the exponents a, (3 governing the 
branching statistics may well depend on the details of the 
random potential, in particular on the precise distribu- 
tion of potential heights and on the two-point correlation 
function of the potential (see Eq. (|l|)). The overall for- 
malism, on the other hand, is more generally applicable 
provided the log-normal behavior of Eq. (0) is satisfied. 
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